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In recent years, a new approach to the theory of nuclear reactions leading to a break-down of 
the interacting subsytems into various channels has been developed. This approach was named 
the Antisymmetrized Molecular Dynamics (AMD), and its main idea lies in the description of 
the nucleons by wave packets in which the antisymmetrization effects (but not other quantum 
effects) are accounted for. In the present review, the basic principles of AMD are illustrated on the 
examples of simplest nuclear systems, and the results are compared with those provided by an exact 
quantum-mechanical description in the Fock-Bargmann space. The applicability region of AMD 
is discussed, in particular, in the cases of systems with discrete spectrum, and a relation between 
the classical AMD trajectories and the quantum distributions is established. At the same time, 
a new interpretation of Brink orbitals and Slater determinants built on them as eigenfunctions of 
the coordinate operator defined in the Fock-Bargmann space is proposed. It is shown that these 
functions form the cluster geometry of nuclear systems. 
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1. INTRODUCTION 

High-energy heavy-ion collisions lead to aJpreak-down 
of nuclei and scattering of nuclear fragmentSLl. This pro- 
cess appears to be essentially non-equilibrated, with a 
large number of nucleons participating in it. Therefore, 
a necessity to apply the kinetic theory to the heavy-i^ 
collision phenomena has been realized some time agi 
However, whereas the standard statistical approach is 
based on classical mechanics, the kinetic theory of nu- 
clear collisions should be based on the quantum equa- 
tions, or rather such a classical limit of these equations 
which would reproduce at least the most important quan- 
tum features of the nuclear dynamics. Thus, a version of 
the kinetic approach, taking into account the quantum 
phenomenon of antisymmetrization and its inflnence on 
the nucleonic classical trajectories was proposedd and was 
named the Antisymmetrized Molecular Dynamic (AMD) . 

LateiO, it became clear that the nucleonic trajectories, 
which are the object of study in AMD, are classical limit 
of wave functions defined in the Hilbert space ofrintcgral 
analytical functions (the Fock-Bargmann space)B. A pe- 
culiar property of these functions is that their variables 
are both coordinates and momenta, i.e. they are defined 
in the phase space. We recall that, usually, transition 
from the coordinate tQ the phase space is done via the 
Wigner density matrixQ, but then one has to deal with the 
alternating signature of this matrix, which contradicts its 
physical sense. On the other hand, the density matrix 
in the Fock-Bargmann space built with the account of 
the Bargmann measure is positively defined at all val- 
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ues of coordinates and momenta, which is its advantage. 
Thus the use of the Fock-Bargmann space in order to re- 
late the quantum and classical statistics is well justified. 
In this space, the classical and quantum mechanics are 
nicely related, including the possibility of the evaluation 
of the classical results in the limiting quasi-classical re- 
gion. This is why we pay a special attention here to basic 
facts concerning the Fock-Bargmann space and illustrate 
them with examples of exactly solvable problems. 

In Section 2, the definition of the Fock-Bargmann 
space is given, and the relation between the wave func- 
tions defined in the coordinate (or momentum) represen- 
tation and their images in the Fock-Bargmann space is 
shown. Section 3 is devoted to the discussion on phase 
trajectories of those systems wave functions of which were 
defined in the previous Section. Influence of the interac- 
tion between the nucleons on their wave functions and 
phase trajectories is studied in Section 4. In Section 5, 
the breathing of wave packets usually used in AMD is 
introduced in the equations of their motion. The conclu- 
sions are presented in the final Section. 



2. FOCK-BARGMANN SPACE 

Although the features of the Fock-Bargmann space 
and functions defined in it are not widely known, some 
of the constructions which are often used in theoretical 
nuclear physics are directly related to this space. 
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2.1. The Brink-Bloch orbital 



The relation 



Among the simplest examples of the wave functions de- 
fined in the Fock-Bargmann space, one will be specially 
discussed here; it is the Brink-Bloch orbital, 

MR) = TT-'/' <^^Pi-y + ^(i^^) - (2-1) 

which was originally proposed^ for the analysis of the 
cluster structure of light nuclei. Here, r is the three- 
dimensional coordinate vector of a nuclcon, and R is 
a vector parameter, introduced to minimize the energy 
of the nu cleus, after a trial functions is constructed 
from (2.1). The Brink-Bloch orbital is a generalization 
of the s-orbital of a nucleon in the harmonic-oscillator 
(harmonic-oscillator) field. As the units of length, mass 
and action we have chosen the oscillator length, the nu- 
cleon mass and the Planck constant h. 

The complex vector R, in general, consists of the real 
and imaginary parts. 



-(R) 



.(R) 



(2.5) 



is now evident. As it should be, the spectrum of the 
operator f is continuous only, and it is due to this fact 
that all its eigenvectors are normalized to the (5-function. 

In addition to all said above, there is one more point 
outlining the special role of the Brink-Bloch orbital in 
the Fock-Bargmann space. It is also the kernel of the 
integral transformation between a wave function ^'(r) of 
the coordinate space and a corresponding wave function 
<i?(R) of the Fock-Bargmann space. 



$(R) = J 0r(R)«'(r)dr. 



(2.6) 



In many cases the integration (2.6) can be done analyti- 
cally. 



2.2. Overlap integral and density matrix 



R 



V2 ' 



R* 
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V2 ' 



(2.2) 



where i is the imaginary unit, ^ is the radius-vector, and 
fj is the momentum. We use the Brink's notation S = R* 
here. 



The orbital (2.1) is the integral analytical function of 



three complex projections Rx,Ry,Rz of the vector R. 
Therefore, it may be considered as a wave function of 
these variables in the (Hilbert) Fock-Bargmann space. 
Then, two questions are to be answered, "how it is nor- 
malized?" and " isn't it an eigenfunction of an operator in 
this space?" The first question is solved by the following 
definition, 



/(r,r') 



:,(R)0r(R)exp{-(RS)}^ 



where the improper integration is over all phase space, 
and 



exp{-(RS)} 



d^dff 



is the Bargmann measureB. The factor {2tt)^ in the de- 
nominator is due to the number of all quantum states in 
the phase space. It is easy to see that 



/(r,r') ==<5(r-r'). 



(2.3) 



Thus, the functions (/)r(R) with different values of r are 
orthogonal, and the components of the vector r are eigen- 
values of these functions. Th e ex plicit form of the oper- 
ator f follows the expression (2.1), 



r = ^(R + VK) 



(2.4) 



After the integration of the eigenfunctions 

0:(R)</.r(R) 

over all their eigenvalues, we arrive to an important ex- 
pression, 

y"<?!);(R)<?!>r(R)dr = exp(RS), (2.7) 

which is often called "the overlap integral" and which, 
also, serves as the density matrix in the Fock-Bargmann 
space. Here we mean the density matrix of pure states, 
and we want to convince ourselves that the expression 
( p.TD satisfy all requirements to be a density matrix. 

The diagonal elements of the standard single-particle 
density matrix defined in the coordinate space (when 
its two vector arguments are considered identical) 
yield expressions for the probabili ty d istribution of the 
coordinatesO. The density matrix (2.7), too, defines the 
probablility distribution, but in the phase space {^, ff} 
and not before being multiplied by the Bargmann mea- 
sure. The resultant expression for the probability distri- 
bution is 



exp(RS)-exp{-(RS)} 



d^df] d^df] 
(2^ " (2^ 



(2.8) 



and we arrive to the uniform distribution over the coor- 
dinates and momenta with the constant density equal to 
unity. Later on, in the discussion on the multi-particle 
wave functions we shall show how this distribution is af- 
fected by the antisymmetrization. 

Of interest is the probability distribution for the state 



:(R)0r(R) exp{-(RS)} - 7r3/2 exp{-(r - ^f} 



(2.9) 
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This distribution is uniform over the momenta and Gaus- 
sian over the coordinates, centered at r. When the orbital 
is used as a trial function in variational calculations, 
^ becomes the variational parameter, Eq. ( |2.9|) is then 
treated as the probability distribution over r, and the 
variational calculation is reduced to a search for an opti- 
mal location of the center ^ of the Gaussian distribution. 
In the classical limit, when Irl ^ 1, i.e. Irl is much 



greater than the oscillator length, the distribution (2.£) 
is reduced to (5(r — f), as it should be for a distribution 
in the state wit h a given value of the radius- vector. 

The relation (2/7) between the density matrix and its 
expansion over a basis of-prthonormalized states illus- 
trates the well-known facttj of diagonality of the density 
matrix in the energy representation. The expansion of 
the density matrix over the eigenstates of a Hamiltonian 
remains diagonal for any Hamiltonian. 

For instance, the density matrix can be expanded over 
the eigenstates of the Hamiltonian of free motion, i.e. 
over the eigenstates of the momentum operator 



k = --(R-V„) 



(2.10) 



The orthonormalized eigenfunctions of this operator, i.e. 
plane waves, labeled by the momentum fc, take the form 



0k(R) =^-'/''exp{ 



R2 



i%/2(Rk) + — }, (2.11) 



which can be proved by solving a first-order differential 
equation with the explicit form of the operator k. An- 
other way to get 0k(R-) is to perform the integration (2.(; ) 
with the plane wave 

*(r) = (27r)3/2exp{-i(kr)}, 

defined in the coordinate space. 

(/.k(R) = y"0r(R.)(27r)3/2exp{-i(kr)}dr. (2.12) 



The expansion of the overlap integral over the states 
with definite values of momentum is similar to (2.7), 



exp(RS) 



UR)0k(R)dk, 



(2.13) 



and the probability distribution for the plane wave (2.11) 
is 

</)^(R)</.k(R)exp{-(RS)} = 7r3/2exp{-(k- ^)2}. 

(2.14) 

Again, this is a Gaussian dependence, but this time for 
the momentum. In the phase plane, the distribution 
function reaches its maxima on the line rf = k, i.e. on 
the classical phase trajectory of the free moving particle 
with momentum k. 



It is worth saying that the Wigner density matrix in 
the same state with the energy E — k^/2 takes the form 



Pw{£.,'n) = (5(k-?7). 



(2.15) 



Therefore, it corresponds to the classical limit of the 
Fock-Bargmann density matrix. 

Let us now turn to the example of the harmonic os- 
cillator which has a discrete spectrum only. In order to 
simplify the situation while keeping valid all its main fea- 
tures, we consider the one-dimensional case first. Then 
R and S are scalar complex variables. 



R 



1 



a/2 



(e + ^77), S^R* 



1 



V2 



where ^ and rj are the coordinate and momentum, re- 
spectively. The general one-dimensional density matrix 
exp(i?5'), similarly to (2.7), can be expanded over the 
states with either given coordinate x, or given momen- 
tum fc, or over the one-dimensional harmonic-oscillator 
states. We are interested in the latter expansion, 



exp(i?S') 



OO 



1 



(i?5)". 



(2.16) 



Eq. (2.16) is followed by an expression for the orthonor- 
malized (with the Bargmann measure) wave fupctions 
(j>n{R) of the one-dimensional harmonic oscillatorllil, 



MR) 



R'\ 



(2.17) 



where n is the number of excitation quanta. As for 
the harmonic-oscillator Hamiltonian -ffosc in the Fock- 
Bargmann space, it is 



(2.18) 



Due to Eq. (2^), the coordinate operator is 
1 _ d 



V2 



(R 



dR 



(2.19) 



and due to (2.10), the one-dimensional momentum is 

d 



and the kinetic energy operator is 



- 1~, 1, 
2 4^ 



dR' 



(2.20) 



(2.21) 



Evidently, 



Hn 



T 



The probability distribution for the harmonic- 
oscillator with the number of quanta n takes the form 



nl 



rjR)MR) eM-RS) 



exp{- 



-}• 



(2.22) 
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It takes only positive values and is concentrated around 
the circle 



and D^'^ are Wigner D-functions. Besides, 



(2.23) 



which is the phase trajectory of the classical oscillator 
with the energy n. Evidently, the left-hand-side part of 
Eq.(2.23) is the classical Hamilton function of the one- 
dimensional oscillator. 

Let us consider in more detail the limiting case of large 
n. At n ^ 1, the following asymptotic formula holds, 



1 



'27rn 



exp ■ 



1 

2n 



.(2.24) 



This is the quantum distribution function for the one- 
dimensional harmonic-oscillator at n ^ 1. It makes us 
possible to calculate the average value of the classical 
Hamilton function and its dispersion. 



^2 + jj2 



= n, 




-n = 



In order to underline the problems with the Wigner 
function Pw{^, »?), we define it following the usual way for 
the harmonic-oscillator with n = 1. Then, 



Pw(tv) = ^ exp{ 



}■ (2.25) 



The latter distribution takes negative values inside the 



circle ^ 



1, which, of course, is a drawback. At 



large values of + though, the distributions pi{£,, if) 
and Pw{^,v) s-re identical. 



2.3. Three-dimensional oscillator 



In the three-dimesnional harmonic-oscillator case, the 



distribution ( 2.16 ) is somewhat modified, 

exp(RS) = ^ — (RS)". (2.26) 

Tl = 

Again, n is the number of excitation quanta. However, 
now (RS)" are wave packets of basis states with the in- 
dices of SU(3) symmetry (n, 0). In order to unambigu- 
iously identify these basis states, we introduce two more 
quantum numbers, the angular momentum / and its pro- 
jection m. Then, 

oo -. 

exp(RS) ^Y.T. Wr.j-{RSrD',^^{n)Dl,^o{n*), 

(2.27) 



n=0 L\ 



where the diagonal elements of the density matrix 

n\{2l + l) 



{n ~ l)\\{n + I + 



(2.28) 



+ Ry + R^ 



R — {Rx, Ry, Rz)- 



Therefore, the product i?"Z?Q„(ri) is the homogeneous 
harmonic polynomial of the nth degree. 

We note that the following condition is satisfied, 

Wr^^l = 1, 

as long as the parities of I and n are the same and I < n. 
In the limiting case of n ^ Z ^ 1, 



21 + 1 , + 

exp{ }. 

n 2n 



(2.29) 



We have arrived to the density matrix for the states 
with n quanta and the angular momentum /. It has the 
form similar to the Gibbs distribution for the rotator with 
the momentum of inertia n and kT = huj, where lu is the 
harmonic-oscillator frequency. 



2.4. Two three-dimensional oscillators 

Two three-dimensional oscillators produce a manyfold 
of such diagonal matrix elements of the density matrix, 
which is better classified by the reduction SU(3) x SU(3) 
^ SU(3). 



exp{(RiSi) + (R2S2)} 



00 00 ^ 

E E ^^(RiSi)"HR2S2)"^ (2.30) 



ni— n2— 



ni\n2 



Then, the product of wave packets with SU(3) symmetry 
(ni,0) and (n2,0) is convenient to present as a super- 
position of normalized wave packets Flnit)(ri^'^o)''^^ ^ with 
the SU(3) symmetry (rii + n2 — 2/i, /i). 



ni!n2! 



(RiSi)"i(R2S2)^ 



E 



'^''^(ni,0)(n2,0) '^(ni,0)(n2,0) ' 



(2.31) 



ni!n2!(ni + ^2 - 2^ + 1)! 

fi\{n - fiy.{n2 - pY-ini + rj2 - /i + 1)! ' 



(2.32) 



where N{nin2) is determined by the normalization con- 
dition, 



min(ni ,712) 

ETT/("l+"2-2/i,A') 
^^(«i,0)(n2,0) 

H=0 



= 1. 



(2.33) 





If ni, 712 ^ A* ^ 1, the expression for the diagonal matrix 
elements of the density matrix is essentially simplified, 



(ni+n2-2At,/j) 



(ni,0)(n2,0) 



1 / nin2 



exp 



nin2 
111 + n2 



(2.34) 



Again, like in the previous subsection, we have arrived 
to a conventional in statistical physics form of the dis- 
tribution function, but this time, for the states with var- 
ious SU(3) symmetry. It is equally important that the 
density matrix elements are squared Clebsch-Gordan co- 
efficients of the SU(3) group, whilst the limiting relation 
( 2.34 ) provides the asymptotic form of these coefficients. 



2.5. Antisymmetrization effects 

Keeping the story within the limits of pictorial one- 
dimensional systems, let us turn to the influence of the 
antisymmetrization to the probability distribution. The 
simplest case woul d be two identical particles. With the 
use of two orbitals (2.1), we perform the antisymmetriza- 



tion and separate out the center-of-mass motion, thus 
getting the following orbital in the c.o.m. frame. 



1 



^/2 



[MR)-M-R)): 



(2.35) 




FIG. 1: Probability distribution p^{^,i]) for the singlet two- 
nucleon system, at (a) k — 0.5, (b) k = 1.0, (c) k = 2.0, 
and (d) k — 3.0. A brighter area corresponds to a highter 
probability. 



_ Xi - X2 _ Ri - R2 

^/2 ' ^/2 ■ 

" {R) is the eigenfunction of the operator 



(2.36) 



with the eigenvalue . The probability distribution for 
this antisymmetric wave function in the Fock-Bargmann 
space is 



r.*{R)rAR)^M~RS) 



= tt"^/^ (cosh(2a;^) ~ cos(2a;r;)) exp(-x^ - 



(2.37) 



Similar transformations occur due to the antisym- 
metrization in the plane wave case, 

0^(i?) = -^(0fc(i?)-<^fc(~i?)), (2.38) 

as well as in the free motion case, 

Pkii^v) - </>fc*(i?)0fe (i?)exp(-i?5) (2.39) 




b) 




d)" 



-2 2 4 




-4 -2 2 4 



-4 -2 2 4 



= 7r-i/2 (cosh(2A:77) - cos(2A:^)) exp(-fc2 - rf). 

Only at fc ^ 1 does the latter distribution limit to a 
superposition of the previously found distribution Pk{R) 



FIG. 2: Probability distribution p^{£,,ri) for the triplet two- 
nucleon system, at (a) k = 0.5, (b) k = 1.0, (c) k — 2.0, and 
(d) k = 3.0. 
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FIG. 3; Density matrices for the singlet (a) and triplet (b) 
two-nucleon systems (see Eqs. (2.41-2.42) for definitions). 



which favours small values of ^,77, and has the same 
asymptotic limit of 1/2 (Fig. 4(b)). 

The probability distribution for the antisymmetric 
harmonic-oscillator states with the number of excitation 
quanta 2n + 1 is 



P2n+l(C,?y) 



1 



(2n+l)!V^^j '"P^ 



(2.43) 



This function reaches its maxima incide a ring, which at 
large n shrinks to the circle 



= 2n + 1, 



which is the classical trajectory at the energy 2n + 1. 



and p-k{R)- At low values of k an ry-dependence appear. 
The surfaces (Ci v) for several values of k are shown in 
Fig. 1. 

The symmetrization of the wave function of two iden- 
tical particles leads to the following expression for the 
probability distribution, 

p+(^, ry) = 7r"i/2 (cosh(2fc?7) -I- cos(2fc^)) exp(-fc2 - rf). 

(2.40) 



The surface described by Eq. (2.40) is shown in Fig. 2. 

Thus, both the symmetrization and the antisym- 
metrization significantly affects the behavior of the dis- 
tribution functions in the vicinity of the origin at low 
values of the free motion energy fc^/2. In particular, if 
ry = 0, and |^| increases, both functions oscillate around 
the value 27r~^/^ exp(— fc^) with the period 27r/fc, remain- 
ing positive. 

The density matrix of the antisymmetric states multi- 
plied by the Bargmann measure is the result of the inte- 
gration of p^{^,ri) over fc, or p^{^,r]) over x, 

w+{^,i]) = J p-{£.,'n)dx = smh{RS)exp{~-RS) 



1 - exp(-2i?S') 



(2.41) 



The antisymmetrization suppresses the probability dis- 
tribution in the vicinity of the origin, but in the outer 
region the probability limits to 1/2 (Fig. 3(a)). 

The integration of pJ(C, ^7) yields the density matrix of 
the symmetric states multiplied by the Bargmann mea- 
sure. 



W+(C,77) 



cosh{RS) exp(-i?S') 



1 + exp(-2i?S') 



(2.42) 



2.6. Three one-dimensional fermions 

Consider now a system of three one-dimensional 
fermions with parallel spins. Let xi^X2, X3 be their coor- 
dinates. Besides, let 



1 2 f X1+ X2 

yi = "^(^1 " ^2), y2 = ^-^yx3 — 



We replace the three independent Fock-Bargmann vari- 
ables i?2, R3 by 



P = \l ^{Ri - R2) = Acosj ^ + i77i, (2.44) 



2 / i?i + i?2 , . . ^ , . 

3 V 2 ' " ^^'"^ = 6 +«??2, 



(2.45) 



Then, the overlap integral of two Slater determinants, 

K,y2{P,Q) and (P, Q), 

constructed of Brink-Bloch orbitals in the c.o.m. system, 
takes the form 



dyi / dy2(t)*y^^y^iP,Q)(t)y,^y^iP,Q) 



j2 (AA*)"2 ■ cos(6m + 3)7Cos(6m -I- 3)7* 

n m— 



(n — 6m — 3)!!(n + 6m + 3)!! 



+ 



^ 'y^' (AA*)"2 • sin6m7sin6m7* 
(n - 6m)!!(n -I- 6m)!! 

n m—l 



(2.46) 
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The expansion (2.46) yields the diagonal elements of 
the density matrix in the harmonic-oscillator representa- 
tion. 



^ ^ {AA*r cos(6m + 3)7 cos(6m + 3)7* ^2 47) 



(ri-6TO-3)!!(n + 6m + 3)!! 



{AA* )" sin 6m7 sin 67717* 
(n — 6m)!!(n -|- 6to)!! 



(2.48) 



The number of quanta n takes only odd values begin- 
ning from n = 3 due to the Pauli principle. Due to the 
same principle, the factors at 7 (7*) in the arguments 
of sines and cosines are proportional to three. Again, if 
ft ^ m 3> 1, we arrive to a simple limiting form. 



dyi / dy2(t>*y^^y^{P,Q)(t>y^^yAP^Q) 



sr^sr^ (AA*)" 2 



(6m -1-3)^ 
2n 



+ 9(PP*)(QQ*)^ - 9([PQ][(P*Q*])(QQ*)}. (2.51) 

These terms are related to the overlap integrals of two 
SU(3) irreducible representations with the indices (0,1) 
and (3,0), respectively. It is worth to note that the (1,1) 
irrep is absent. 

The first terms of the expansion of the overlap inte- 
gral determine its behavior at small values of the vectors 
P, Q; P*, Q*. Thus, the very first term, 

i([PQ][P*Q*]), 

is the overlap of the translation-invariant wave functions 
of the Elliott's SU(3) modelll3 of three fermions with the 
same direction of their spins, while the fourth power of 
this term is the overlap for ^^C, in the same model. If 
the vectors P and Q are collinear and arc of the same 
length, we can set P = Q = R. Then, the first term of 
the expansion (2.51) vanishes, and the second ttakes the 
form 



(RR*)^ 



which is responsible for the linear structure of the wave 
function. Again, this term corresponds to the overlap in- 
tegral of the three spin-correlated fermions but occupying 
the states 



X cos(6to -I- 3)7Cos(6m -|- 3)7* 



EE 



{AA*y 



■ — cxp 

n 



2n 



(2.49) 



sin 67717 sin 6m7* 



[0,0,0] [1,0,0] [2, 0,0] 

of the s-,p- and sd— shells, respectively. Here we use 
the traditional notation [ti^,, Uy, Uz] for the single-particle 
harmonic-oscillator states with the numbers of quanta 
along the axis i. 



which is similar to the one for the three-dimensional os- 
cillator. 

With the density matrix wl^ ^i^i^ V1j^2,V2) at our dis- 
posal, it is easy to find the single-particle density matrix 
mCCii'/i) the form of the integration. 



(2.50) 



w'n,m(6,'7i) = / <,™(6,?7i;6,?72) 



X exp{- 



}- 



2tt 



2.7. Same system, three-dimensional case 

In the case of three three-dimensional fermions having 
the same direction of spin, the first two terms of the 
expansion of the overlap integral take the form 

i([PQ][P*Q*])+ 



24 



{(PP*)^ - 3(PP*)(QP*)^ - 3(PP*)(PQ*) 



PHASE TRAJECTORIES 



The orbitals (2.1) clear the way of determining the 
classical phase trajectories of the wave packets of a sys- 
tem of A fermions. In practice, in order to derive the 
equations of classical dynamics and find the phase trajec- 
tories, several stages are needed to be gone through. At 
first, the Slater determinants are constructed of the or- 



bitals (2.1) thus meeting the Pauli principle reqirements. 
Then, these determinants are declared as the trial func- 
tions of the variational problem, with the complex vectors 

R,, z = l,2,...A 

being the variational parameters depending on time t. 
Next stage is the calculation of the Lagrange function. 
Finally, the least action principle is applied, which results 
in the classical equations for 



Kit) 



V2 



and, in the end, the phase trajectories are found. 
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At a given energy, each trajectory of the system of 
fermions is located on a hyper-surface in the phase 
space. It is completely defined by the integrals of mo- 
tion. On the other hand, in the Fock-Bargmann space, 
the squared absolute value of the wave function of the 
system multiplied by the Bargmann measure gives the 
probability density distribution in the same phase space. 
By comparing the phase trajectories and the quantum 
distributions, we may judge, to what degree and under 
what conditions the phase trajectories reproduce the real 
situation. 



3.1. Motion in one dimension 

It is especially easy to find the cl assi c trajectories of 
motion in one dimension using Eq. (2.1). The following 
relation is used, 



KmmAR)) (3.1) 



{<j,l{R)\UR)) 



Here H is the quantum Hamiltonian of a one-dimensional 
system in the coordinate representation, 'H{(,,ri) is the 
Hamilton function defined as the result of variation of 
the Hamiltonian over the Brink-Bloch orbitals. Thus, for 
the one-dimensional free motion, when the Hamiltonian 
is 



H 



the Hamilton function appears to be 



(3.2) 



It follows from (|3.2D that the wave packet (^T|) moves 
as a free particle having momentum rj and the mass of 
the nucleon. However, as E — T-C{^,ri), the total energy 
of the wave packet contains an additional term equal to 
— 1/4. It is this energy which is required to create a one- 
dimensional wave packet. As any wave packet, this one 
cannot exist for a long time and should spread apart. As 
far as we neglect this phenomenon and do not consider 
its mechanism, the minimal energy of the wave packet 
appears to be 1/4. It remains so even when an attractive 
potential is switched on and a spectrum of states with 
positive energy not more than 1/4 appears. Naturally, 
the quantum description does not have such a drawback. 

As for the phase trajectories of free motion (3.2), the 
dispersion relation 



V 



V^E - 1/2, 



somewhat differs from both the classical one and the one 
f ollow ing the line of maxima of the distribution function 

For the one-dimensional oscillator 



H 



'2da;2 



1 



(3.3) 



therefore. 



(3.4) 



and the phase trajectories at energy E satisfy the relation 



2^ ' 



- = E. 
2 



(3.5) 



As expected, at any E the motion is finite while the spec- 
trum is continuous but hmited (E > 1/2). At E = 1/2 
the phase trajectory degenerates into the point of origin 
(^ = 0, 77 = 0). We remind tha t the quantum distribu- 
tion in the same state (see Eq. ( p.22 )) is 



It has a maximum at the origin point. 



3.2. Free motion of triplet pair of fermions 

In order to understand how the anstisymmetrization 
affects the phase trajectories of free one-dimensional mo- 
tion of two fermions in the triplet state, we define the 
classical Hamilton function. 



(3.6) 



where, again. 



H 



2 dx'^ ' 



Simple calculations lead to the following rcsultc 



2 



e 



(3.7) 



Th e sum of two last terms in the right-hand side of 
( |3.7|) is shown in Fig. 4. Having set H^{£,,ri) — E, we 
obtain the phase trajectory at energy E. Note that the 
minimal energy is equal to 1/4. Besides, as 



coth 



1>0, 



the Hamilton function contains a repulsion. Therefore, 
as long as 1/4 < < 3/4, the trajectory begins at large 
positive values of ^ and 77 — yj2E — 1/2 and after the 
refiection does not enter the region of negative ^. There 
appears a turning point where the momentum is zero. 
Having passed through this point, the trajectory remains 
in the region of positive ^ with the sign of momentum 
alternated (Fig. 5). If £' > 3/4, the motion never stops 
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FIG. 4: Effective potential in the triplet fermion pair case. 
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FIG. 5: Phase trajectories in the triplet fermion pair case, a) 
E = 1.5, h) E = 1.0, c) E = 0.7, d) E = 0.4. 
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FIG. 6: Effective potential in the singlet fermion pair case. 



3.3. Singlet pair 

In the case of the singlet pair of free fermions, the 
Hamilton function takes the form 



e 



^ftanh^ 



(3.8) 



The last term in the r.h.s. of ( p^ ) appears, again, due 
to the symmetrization but, unlike in (3.7), it corresponds 
to the attraction in the vicinity of the origin of the phase 
space, as 



tanh 



1 < 0. 



but slows down. Then the trajectory enters the region 
of negative ^ where it can be reproduced by a reflection 
around the ry axis. 

The repulsion visible in the phase trajectories of the 
triplet pair of fermions at small values of ^ and 77 agrees 
well with the conventional wisdom: the Pauli principle 
forbids the fermions having parallel spins to be located 
at the same point. 



The probability distribution (2.40) for the triplet pair 
with energy E = limits to the phase trajectory 

H~(^, rf) = E only ii E ^ 3/4, i.e., when the influence of 
the Pauli principle and the antisymmetrization is negligi- 
ble. It is then when the classical result is the asymptotic 
limit of the quantum one. However, at small energies, 
where the corrections to the classical results due to the 
antisymmetrization become significant, there is no anal- 
ogy between the quantum probability distribution and 
the phase trajectories. In order for this analogy to ap- 
pear, some additional corrections of classical equations 
are required. These corrections are mostly due to the 
spreading of wave packets described by classical equa- 
tions. 



The behavior of two last terms in the right-hand side 
of (3.8) is shown in Fig. 6. The minimum of this sum 
is located at the circle centered around the origin. Hav- 
ing added rf /2 to this sum, we get two minima of the 
Hamilton function at 

e-±eo, Co = 1.13, r/ = 0, 

when Erain = 0.11. In these minima, the phase trajec- 
tories degenerate into a point. Then, when the energy 
exceeds its minimal value, there appear closed phase tra- 
jectories of finite motion (Fig. 7). The motion becomes 
infinite \i E > 1/4. However, an explicit relation between 
the infinit e traj ectories and the quantum probability dis- 
tribution (2.4C) becomes evident only if _E 3> 1/4 and the 
influence of the symmetrization in both the phase trajec- 
tories and probability distributions becomes negligible. 
As for the finite trajectories, they are the consequence of 
the approximation applied at the derivation of the classi- 
cal equations. The approximation under question is that 
the wave packets are considered as having a fixed width, 
although in fact they spread apart as they move. This 
phenomenon may be neglected at high energies only. 
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FIG. 7: Phase trajectories in the singlet fermion pair case, a) 
E = 0.40, h) E = 0.30, c) E ^ 0.25, A) E ^ 0.25, and e) 
E = 0.15. 



4. INTERACTING FERMIONS 

We now introduce a fermion-fermion interaction poten- 
tial in order to study those problems which appear when 
a system has both discrete and continuous spectra. The 
interaction will be simulated by a Gaussian interactive 
potential, 



U{x) — — Vb exp{— ax^}. 



(4.1) 



First, wc find the spectrum of states defined in the Fock- 
Bargmann space, then we perform a transformation to 
the classical equations and determine the phase trajecto- 
ries. 



where z -'^ = 1 + a. A similar overlap with the kinetic 
energy operator has been calculated earlier, 

{S\f\R) - (-\{R - Sf + -\ eMRS). (4.5) 



The overlap integrals (4.4) and (4.E) are the matrix 
elements of the operators involved between the orbitals 
(4.2) and ( |4.3| ). But these orbitals are superpositions 
of all orthonormalized basis functions of one-dimensional 
harmonic oscillator defined in the conventional coordi- 
nate space. The coefficients of these superpositions are 
the basis functions of the same oscillator but in the Fock- 



Bargmann space. For example, for (4.2) 



1 1 

MR) = E ^ gn(x) exp{-xV2}^fl", 

(4.6) 

where Hn{x) are Hermite polynomials. Having inte- 
grated over X, we found a linear combination of all matrix 
elements. It is our task now to extract the proper matrix 
elements from this combination. 

A matrix element {S\0\R) of an operator O can be 
expanded as follows, 

{S\d\R) = J2T. -F=R''{n\^\n)^S^ 



n=0 h=0 



(4.7) 



where 



(n\6\n) = 



dx 



4.1. Solution of the quantum wave equation 

In the framework of quantum approach, we find the 
overlap integral 



{S\U{x)\R) = / MS)U{x)MR)dx 



of the potential energy operator (4.1) between the one- 
dimensional Brink-Bloch orbitals 

</.,(i?)=^-i/4exp{-y + V2i?x-^}, (4.2) 



/),(5)=^-i/'*exp{-y + x/25x-^}. (4.3) 



It is easy to see that 



z- 1 



{S\U{x)\R) = -z^/^Voexp{zRS + ^-{R^ + S^)}, 



(4.4) 



xi/„(x)exp{-xV2}0- 



1 



,H.n{x)exp{-x^/2}, 



which shows the way to solve our task. Indeed, as 

oc oo -. -J 

{S\U{x)\R) -EE -^R''{n\U{x)\fi)^S\ 



(4.8) 



the partial matrix elements are 

inin(n, 

{n\Uix)\n)^~Vo E 



m=0 



m\{n — my.\(h — m)!! 



(4.9) 



Eqs.(4.5) and (4.7) are followed by three well-known 
expressions for the matrix elements of the kinetic energy 
operator of one-dimensional oscillator, 

1 



(n + 2\T\n) = --^{n + l){n + 2), (4.10) 
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{n~2\T\n)^--^ 



(nlTln) = —n -i — . 
\ I I / 2 4 



l)n, 



(4.11) 



(4.12) 



We used the matrix elements between the Brink-Bloch 
orbits and basis functions defined in the Fock-Bargmann 
space. We have specificaUy considered a simplest case 
when all calculations are free of background details. How- 
ever, the algorythm outlined above remains applicable in 
a more complicated situation if the overlap integrals and 
the basis functions in the Fock-Bargmann representation 
arc known. 

Now, using the matrix elements of the Hamiltonian 
found above, we consider a set of algebraic equations 

oo 

Y,{{n\H\n} - M„3)C„ = 0, n = 0, 1, 2, 

h 

(4.13) 

for the coefficients Cn of the expansion of the wave func- 
tion 'I'(-R) over the harmonic-oscillator basis. Naturally, 



We can analyze an example of free motion, when 
U{x) = 0. As in the general case, the continuum is de- 
generate in the two-fold way, and in order to remove this 
degeneracy, parity is introduced as an integral of motion 
along with energy E. Evidently, the wave function of an 
even state 



(4.16) 



Tl = 



contains only the basis states proportional to even powers 
of R, and the wave function of an odd state 



vI/-(i?) = ^C- 



1 



r!=0 



v/(2" + l)! 



-.R 



2n+l 



(4.17) 



contains those proportional to odd powers of R only. 

The values of the coefficients C+ (C~ ) of the expansion 
of the wave function of free motion follow the expression 



^^(^R) ^ Tr~'/Uxp{-^ - iV2Rk + ^} (4.18) 



^{R)^Y.Cn^R-. 



(4.14) 



The same coefficients C„ enter the expansion 



which is a one-dimensional analogue of the plan e wave 
( ^.11 ) with the momentum k. Note that ( 4.18 ) is not 
only the wave function of free motion in one dimension 
but also the generating function of the Hermite poly- 
nomials ii™(fc), with R being the generating parameter. 
Thereforell^, 



^{x)^y^Cn , , ^ g„(x)exp{-xV2} (4.15) ^k{R)^\^ 
■H. ^/2"r^!^/7^ "=o 



of the wave fun ction in the coordinate representation 
which has ( 4.14 ) as its image in the Fock-Bargmann 

space. 

The solution of the set ( [4.13| ) yields the spectrum of the 
pote ntial U (x). Sometimes the applicability of the expan- 
sion ( 4.14 ) to the contimmm states is questioned. This is 
related to the fact that, in the coordinate representation, 
one-dimensional continuum wave functions oscillate with 
X increases with a non-decreasing amplitude. Meanwhile, 
the harmonic-oscillator basis functions fall sharply with 
the increase of cc, and the question appears, whether it is 
possible to expand a non-decreasing oscillating function 
in a series over the harmonic-oscillator functions. Pre- 
cisely speaking, the convergence of such a series is under 
question. 

In the defen ce of the methode, we make a note on 
the expansion (p^ , which is a Tailor series of an en- 
tire function ^'(-R). By definition of a whole function, 
its expansion over powers of R coverges everywhere in 
the complex plane except the infinity where this func- 
tion has an irregularity. This general statement on the 
convergence of ( 4.14 ) is confirmed by the behavior of the 



;i/„(fc)exp{-^}^i?». 



c+(fc) 



(-ir 



,/22»(2r 



= iJ2n(fc)exp{ -}, 

'TT ^ 



(4.19) 



(4.20) 



C-{k) = , ^ ' ^g2n+i(fc).exp{- — }. 

V22"+i(2n + l)!V^ 2^ 

(4.21) 

We underline some important facts here. First, as ex - 
pected, the coefficients C+ and C~ satisfy the aet ( 4.15 ) 
at U{x) = 0. This is a consequence of the knowEll3 recur- 
rent relations for the Hermite polynomials which are just 
identical to the set (4.13) in this particular case. Second, 



these coefficients are the eigenfunctions of the harmonic 
oscillator defined in the momenturp-representation. 

Finally, the asymptotic behavioJij of C+ (and C~ ) in 
the limiting case n ^ 1 



coefficients of the series: |C„| are limited and l/Vn! con- 
verges to zero. In a circle of any radius centerd around 
R = Q this convergence is uniform. It ceases to be steady 
out of this circle. 



Cni^) - (-l)"^n-l/4c0s(fcV4^ 



1), (4.22) 



'i/''sin(A:V4n + 3) (4.23) 
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is identical to that of the solution of the one-dimensional 
Schrodinger equation in the coordinate representation for 
free motion with energy E — Here, the role of the 

absolute value of the coordinate is played by the value 
of the turning point of motion in the harmonic-oscillator 
field, i.e. y/'in + 1, if energy is 2n + 1/2, and y/An + 1, 
when it equals 2n + 3/2. The expansion coefficients are 
normalized. 



(4.24) 



Naturally, the potential U{x) affects the values of the 
expansion coefficients but we still can define a priori its 
asymptotic behavior which takes a familiar form. 



C+(fc) - (-l)"^n-i/-* cos(fcV4n + 1 + S+{k)), 



(4.25) 



Vq is approaching zero, then — e — !■ 0, and the wave func- 
tion becomes slowly falling and delocalized, spreading far 
away from the origin. 

The density probability distribution Pt{^,r]) in the 
phase space is defined in terms of a two-fold summation. 



n=On=0 



X Q 



}■ 



(4.28) 



As u sua l and as opposite to the classical result ((3.7) 
and ( |3.8| )) for the phase trajectories, the infinite motion 
is possible at any E > 0. 



4.3. Phase trajectories for the potential (|4.l]) 



Cnik) - (-l)"^n-^/^sin(/^:^/4^^^3 + <5-(fc)) 

(4.26) 

Thus, at any n > hq, when the asymptotic expressions 
are defined with a pre-defined precision, all coefficients 
are defined by one unknown parameter only, the phase 
shift S^{k) {S^{k)), which has to be found along with 
the coefficients (Cr7)i at < no. In this way, the 
infinite set of equations is reduced to a closed set of no + 
1 equations. The choise of no depends on the desired 
precision of calculations. 

These general considerations basic speculations remain 
valid in the three-dimensional case, as well as in the cases 
of many particle systems, and with the Coulomb inter- 
action included. In any case, it is necessary to detemine 
the explicit form of the asymptotic behavior, which, in 
general, contains the scattering K- or 5-matrix elements, 
and express it in terms of quantittn numbers of the muti- 
dimensional harmonic oscillatorliJ. 



4.2. Discussion on solutions 



Let the attractive potential (4J) generates a single 
bound state with energy — e, positive parity and the ex- 
pansion coefficients {C^}. As for the continuum states, 
they are labeled by the momentum at infinity k and par- 
ity. The expansion coefficients are {C+(fc)} for even (sin- 
glet) and {C~{k)} for odd (triplet) states. 

Energy and structure of the wave function of the 
ground state depend on the depth Vo of the potential. 
If Vb > 1 then 



(4.27) 



and the wave function is localized around the origin in 
both coordinate and phase representations. If, however. 



If motion of a particle in the field of the potential 
U{x) = —Voexp{—ax^} is studied, when the Pauli prin- 
ciple does not affect the results in any way, the following 
term should be added to the Hamilton function of free 
motion (^^), 



{rjR)\u{x)\MR)) 

{ct,l{R)\MR)) 



-zVVoCxp{^(i? + 5)n. 



(4.29) 



Then, 



m, ^) = y + i - ^'^'^0 °'^Pi-(i - ^)^'}- (4-30) 

As in the free motion case, infinite trajectories appear 
only if E > 1/4. Note that the shape of the potential 
remain Gaussian after averaging over the Brink-Bloch 
orbitals. However, as 1 > z > 0, the depth of the resul- 
tant potential appears to be ^/z times less than that of 
the original potential while its width increases ^/z times. 
Therefore, even in the limiting case of high energy when 
the term 1/4 in the expression for the Hamilton func- 
tion can be neglected, the scattering picture should dif- 
fer from that provided by the potential (4.1). There- 
fore, when calculating the trajectories of infinite motion 
in AMD, some corrections are required to produce ad- 
equate results. In fact, such corrections are introduced 
in heavy- ion collision studies. As it became clear now, 
a modification of a nucleon-nucleon potential used (usu- 
ally, a Volkov potentialla) cannot yield correct results 
neither at low energies (classical trajectories do not even 
resemble the behavior of the quantum wave function), 
nor at high energies (the potential is deformed). Instead 
of such a modification, an^ditional two-nucleon scatter- 
ing is introduced in AMDEj. Actually, then the nucleon- 
nucleon scattering is taken into account twice; once via 
the Volkov potential and once via the additional scatter- 
ing cross-section adjusted to the experimental value. 



13 



The minimal energy of finite motion of the particle 
corresponds to i? = 0. Here the phase trajectory de- 
generates into the origin point. In the coordinate rep- 
resentation, the AMD wave function is reduced to the 
harmonic-oscillator ground state, 



/.:(i? = 0) = 7r-V4exp{--}. 



(4.31) 



This approximation can be satisfactory if the value of 
the oscillator radius is optimal and the potential is deep 
enough. However, it looses its precision as Vq decreases. 



4.4. Projection in AMD 

The projection onto the states with a definite values of 
the angular momentum presents a problem in AMD. As 
an example, consider the motion of a three-dimensional 
particle in the central potential field. 



C/(r) 



-Vo exp{— ar^}. 



(4.32) 



The standard Hamilton function in the phase space (r, p) 
then has the form, 



(4.33) 



where pr is the radial projection of the momentum, M — 
[rp] is the angular momentum. The projection procedure 
starts from fixing the value of M. Then a point r = vm 
is determined in which the function 



A/f2 

F(r,M2) = — + U{r) 



(4.34) 



reaches its minimum, and the minimal energy at given 



Vq exp{-ar^^} 



(4.35) 



is calculated. There is a maximum of F{r, M^) at rmax! 
the maximum value -F'max(-^^) is positive. 

In the plane spanned by (r, p^), the phase trajec- 
tory shrinks to the point at the energy Eq{M'^). If 
Fma^M^) > EiM"^) > EoiM"^), the phase trajectories 
are closed, concentrated around the point r = tm, Pr — 
0. If E{M^) > Fjna^iM"^), the motion becomes infinite 
and the phase trajectories go to the infinity. 

The Hamilton functi on Ti . of the wave packet in the 
field of the potential (4.32) differ from (4.33) in some 
details only; 

n = ^ + l- z^/Vo exp{-i^en, (4.36) 



These details are the energy 3/4 needed to create the 
wave packet, the factor z'^/^ and the factor 1 — z instead 
of a in the exponent. There are no other differences, 
though. The minimal and, thus, the ground state, energy 
is reached if both t^^ and M are zeros. Then the minimum 
corresponds to the zero value of ^, and that gives 



En = 



which means that a degenerate to a point trajectory ex- 
ists at any small positive Vq. This contradicts the state- 
ment of quantum mechanics that a bound state exists 
only if Vq exceeds some critical value. 

Furthermore, the Hamilton function ( 4.36| ) provides us 
with a continuous spectrum of finite states, and some 
additional conditions are required to select those states 
which do not contradict quantum mechanicaJ princi- 
ples. The BohfySommerfeld quantization ruleaJ or their 
generalization ^ may serve as such conditions. 

Besides, Eq.( 4.36 ) allow for not only integer values of 
the angular momentum M, but any other, too. In order 
to overcome this difficulty, the Brink-Bloch orbital has 
to be projected onto a state with the definite value of the 
angular momentum. This procedure is commonly known 
as projection before variation. Let us trace its realization 
on the simple example of the zero angular momentum 
state. In this case, the overlap integral takes the form 



exp{RSt}dt = 



sinh i?S' 
RS ■ 



(4.37) 



The integration here is over t, the cosine of the angle 
between the vectors R and In fact, we follow the 
known Peierls-Yoccoz methodtS. After the projection is 
done, the Hamilton function is 



Ha 



2 



coth 



1 - exp{-z(^'^ + if)} 
l-exp{-(e2 + r;2)} 



exp 



(4.38) 



where 



R=:^ + ir], S = R* 



The minimization of the Hamilton function (4.38) leads 
to the same ground state energy that we o btain ed ear- 
lier. However, the change of the potential (4.32) by the 
additional repulsive term 

Vi exp{-Q!ir^}, 

quite affects the result. Now, in order to find the ground 
state energy, we need to minimize the function 
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i/2t^ 1 - exp{-z^^} 



exp 



1 - z 



r + 



1/2. , 1 - exp{-zi^^} ^ 1 - ^1 ,^2 



1 _ exp{-C2} 



■ exp 



r (4.39) 



over the semi-axis < ^ < oo. If the projection before 
variation is not performed, the function to be minimized 
is 



^ - z^/Vo exp{-i^en + exp{-i-- 



(4.40) 



In both cases, the minimum is located at different values 
of ^ 7^ 0. Evidently, the projection lowers the minimum 
value. 



4.5. Influence of the Pauli principle to the 
Hamilton function 



Having taken into account the potential (4.1) we 



should modify the Hamilton functions (3.7) and (3 



The additional term for the function of the singlet pair 



{ci,riR)\uix)\4>+m 



i/2,r exp{-(l - z)e} + cxp{-e - V} 

-Z Vq 



1 + exp{— — 7]^} 



(4.41) 



exp{— (1 — z)^^} — exp{— — zri"^} 
1 — exp{— ^2 — ry^} 



(4.43) 



We have already stressed the dependence of the results 
on the depth of the potential in AMD. If the potential 
is deep both approaches (AMD and exact) show similar 
resluts for the ground state energy and the wave func- 
tions (Fig. 8). On the contrary, if the potential is so 
shallow that the exact solution for the ground-state en- 
ergy appears to be just below the threshold, the AMD 
approach yields a positive energy but a bound-like wave 
function (Fig. 9)! The same figure shows what kind of 
function correspond to the ground-state energy of AMD 
in the exact case. 

We would like to note that the results discussed above 
are directly related to the variational calculations per- 
formed in AMD and show the way to improve them. A 
typical object of study for AMD is a light nucleus fea- 
turing an a-cluster structure. This structure is a con- 
sequence of the strong interaction between the nucleons 
constituing the same a-cluster as compared to the weak 
interaction between the nucleons belonging to different 
a-clusters, which is due to the Pauli principle and the 
exchange nature of nuclear force facilitating the satura- 
tion of nuclear matter. We may refer to the well-known 
experimental fact: the binding energy of light nuclei is 
close to the sum of binding energies of their a-clusters. 
In AMD, the binding energy of a-clusters is reproduced 
rather well but the contribution of the cluster-cluster in- 
teraction to the binding energy is underevaluated. 



4.6. AMD and RGM 



Thus, the averaging of the Gaussian potential over the 
symmetrized Brink-Bloch orbitals leads to an interaction 
in the phase space which does not depend not only on the 
coordinate ^ but also on the momentum rj. However, at 
large distances between the particles, when ^ 1, this 
interaction becomes Gaussian, as it was in the previous 
example when the Pauli principle didn't play any role. 

As a result, the Hamiton function of the singlet pair 
takes the form. 



-^---^^-i^^ftanh^ 



H-^(C,.)-y + 4 4 



^i/2y exp{-(l - z)e} + cxp{-e^ - zr?^} 
1 + exp{— ^2 — jy^} 



(4.42) 



In the triplet pair case, the Hamilton functions is 
changed in a similar manner. 



In order to explain the notion of the previos paragraph 
it is worthy to compare AMD to the Resonating Group 
Method (RGM) which also employs Brink-Bloch orbitals. 
Let us consider, for the sake of simplicity, an a-cluster nu- 
cleus consisting of n clusters {A — An nucleons). Whereas 
the AMD wave function of this nucleus is a Slater de- 
terminant of A different Brink-Bloch orbitals, there are 
four times less orbitals in RGM, since four nucleons of 
each cluster are assumed to have identical orbitals. Af- 
ter the c.o.m. wave function is separated and a new set 
of Jacobi-like vector parameters is introduced as a lin- 
ear combination of the original parameters of the Brink- 
Bloch orbitals, AMD provides a method to study both 
nucleonic and cluster motion. Variational AMD calcu- 
lations with the Volkov potential show that the optimal 
values of the vector parameters for nucleons of the same 
cluster are identical while the vector parameters deter- 
mining cluster locations differ from zero. The conclusion 
is made that there is an a-clusterization in light nuclei. 
On the other hand, this justifies the choice of orbitals in 
RGM. 

Surely, this conclusion is valid as long as the Volkov 
force can be considered as describing the real picture. 
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FIG. 8: (Upper panel) Ground state energy of the triplet 
pair of fermions interacting via the deep potential V{x) — 
— Vo exp(— x^), Vo = 83 MeV. The potential is show n by the 



solid thick line. Effective potential (see Eq.(4.43)) for the 
wave packets is shown by the dashed thick line. (Lower panel) 
The wave functions corresponding to these energies. 



FIG. 9: Same as in the previous figure but the potential 
is shallow, Vo = 33 MeV. The short-dashed line shows the 
quantum-mechanical solution at E = 4.58 MeV, belonging to 
the continuum. 



For it is known that this force does not satisfy the satu- 
ration condition, and calculations of sd-shell nuclei with 
it does now exhibit any clusterization even if the cluster 
structure is expected. 

The shape of the nucleon-nucleon function in AMD is 
almost independent on whether the nucleons it describes 
belong to the same cluster or not. This shape is chosen to 



be the best for the a-cluster nucleons which leads to an 
error in the calculation of the cluster-cluster interaction 
energy. In RGM, the wave functions of the cluster-cluster 
motion is defined by an optimization procedure. Thus 
the problems facing AMD are avoided and the results 
are bettered. 

In Ref.0, the RGM orbitals were used in derivation 
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and analysis of the classical equations of motion for the 
problem of scattering of the a-particlc off the ^^C nucleus 
considered as three a-particles. 

It is important to note that both AMD and RGM 
are generalizations of the shell model, or rather, Elliott's 
SU(3) modclB. The latter is the hmiting case in AMD 
and RGM, when all vector parameters limit to zero. 



5. SPREADING OF WAVE PACKETS 

Spreading of wave packets is a well-known phenomenon 
in quantum mechanics. Nevertheless, this phenomenon 
has to be neglected when the problems of the contin- 
uum are studied using dynamical equations of AMD. It 
seems to be justified if the time of interaction of collid- 
ing nucleons is much less than the time of spreading of 
their wave packets. The less the energy spent to create a 
wave packet, the faster it spreads. For a single particle, 
this energy equals 3/4. We use this simple example to 
show how, using classical considerations, to describe the 
spreading of a wave packet and estimate the time scale 
of this process. 



5.1. Wave packet of radial oscillations 

Consider the wave packet 
1 1 



I 1+g 2 

^3/4 (1 ^ £)3/2 <j 2 1 - e ^ 



(5.1) 



of a particle in the three-dimensional space0. It is usu- 
ally used in the description of the radial oscillatory mode. 
The width of the packet depend on the generating param- 
eter e. The wave packet spreads apart if e approaches — 1. 
The overlap integral 



(0,(e)|0,(£)) = 



1 



(l-£e*)3/2 



£_£(^±iZ^,",- (5.2) 



^^^r(3/2)r(n+l) 



provides us with an image 
't>n{e) = 



r(n + 3/2) ^„ 



r(3/2)r(n-Hl) 



(5.3) 



where L^J'^ir^) a re th e Laguerre polynomials. 

The functions ( |5.3D are orthonormalized with the mea- 
sure 



(1 -££*)- 



^il2 dadf3 



, e = a + i/S 



in the unit circle |e| < 1. Indeed 

0:,(e)0„(e)(l 

e|<l 



(5.5) 



Similarly to the orbital {2A), the wave packet (5.1) is 
an eigenfunction of the squared radius operator. 



^2 a 3, 



(5.6) 



defined in the e-representation (inside the circle of unit 
length on the complex plane), i.e.. 



(5.7) 



Later on we shall be considering free motion, with the 
Hamiltonian 



2 



(5.8) 



The same Hamiltonian can be written in the e- 
representation. 



1 , ,T 9 3 , 



(5.9) 



The eigenfunctions of the latter, corresponding to the 
energy E — k'^ /2, are 



0fe(e) = 



1 



1 



7^3/4 (1 + e)3/; 



■ cxp 



1 1-e 



2 1 + e 



(5.10) 



As we had before, we now face a dilemma; either we 
write down the wave equation in the e-representation and 
search for the exact quantum-mechanical solution to the 
problem of the radial motion, or, by studying evolution 
of the wave packet, we reduce the problem to a simpler, 
classical one. Beyond doubt, in the latter case we shall 
get only an approximated description; a phase trajec- 
tory instead of a quantum probability distribution in the 
phase space. Below, we follow this way. 



of the basis functions of the three-dimensional harmonic 
oscillator with the angular momentum L = 0, defined in 
the comple x pl ane of e. That it is so is easy to see having 
expanded (5.1) over the powers of e, 



n=0 



I r(n-i-i) 

2nr{n + 3/2) 



5.2. Phase trajectories of radial mode 

It is easy to see that the classical Hamilton function of 
free motion is 



n 



{Me)\H\Me)) 



3 (l+£)(l + £*) 

4 l-e*e 



(5.4) 



(5.11) 
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FIG. 10; Phase trajectories for the radial mode. Shown are 
the trajectories for x ~ 0-5 = •y/3/4; solid line), x ~ 1-0 
(fc — yJ'i/2\ short-dashed line), and x = 2.0 {k — \/3; solid 
line) . 




FIG. 11: Probability distribution function W{r,t) in the co- 
ordinate space (notmalized to unity at each moment of time). 
The curves correspond to f = (solid), t = 1 (long-dashed), 
and t = 2 (short-dashed). 



Eq.( ^.ll ) is followed by the equation for the phase tra- 
jectory at a positive energy E, 



X 



-)' + /3' = (l- 



X 



(5.12) 



where x — 4i?/3. Each of the trajectories is a circle cen- 
tered around the point —1/(1 + x) located on the real 
axis a (Fig. 10). The radius of the circle is 1— 1/(1 + x). 

As a phase trajectory approach the point a — —1, P — 
0, the width of the wave packet increase, and the packet 
completely spreads apart at this point. As mentioned 
above, the energy spent to create the packet is E = 3/4. 
In this case, x — the radius of the phase trajectory is 
1/2, and it joins the point e = in the complex plane 
with the point a = —1, /3 = 0. 

In order to estimate the time of spreading we turn 
to the classical equations of dynamics. These equations 



follow the least action, principle and the expression for 
the Lagrange functiono. 



(5.13) 



As a result of the variation over e* we arrive to the fol- 
lowing ordinary first-order differential equation for e. 



(5.14) 



or, alternatively, to the set of equations for a and /3, 

d = -(l + a)/3, p=l{{l + af-/3^}. (5.15) 



Integration of Eq.(5.14) with the initial condition e(0) — 
yields 



2t 



(5.16) 



The function (3{t) reaches its minimum at t = 2, when the 
phase trajectory has elapsed the first quarter of the cir- 
cle. It takes infinite amount of time to elapse the second 
quarter. 



Having used Eqs. (5. If), (5.1), and taking into account 
the normalization factor, we obtain the time dependence 
of the probability density of values of r in the coordinate 
representation. 



W{r, t) 



(f)r{e)(t>ris) 

1 1 



7r3/2 (l+t2)3/; 



■ exp 



1 + 



.(5.17) 



Thus, the time dependence of a wave p acket at rest is ex- 
plicitly found. It is illustrated by Eq. ( 5.17 ) and Fig. 11. 

Finally, in order to understand the relation between 
the quantum result in the Fock-Bargmann space and the 
classical result, it is useful to compare the phase trajec- 



tory (5.12) with the quantum distribution in the phase 
plane for the state 0A:(e) with the kinetic energy k'^/2. 



-1/2 ^-3/2 



1 



X exp{- 



1 - ££* 



(l + e)(l + e*) 



3/2 



(5.18) 



{l + e){l + e*) 

5.3. Spreading of Brink— Bloch wave packet 

Finally we show how to modify the expression for the 
Brink-Bloch wave packet in order to completely take into 
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account its evolution due to its spreading as it moves in 
the three-dimensional space. Assume 



-(e,R) 



1 1 

^ (1 - e)3/2 



r 1 l + £ 2 

exp-^ — r 

2 l-£ 



-^(Rr) - }. 5.19) 



The overlap integral of ( 5.1£ ) is 



-(£,R)|,/)r(£,R)) 



(1 - ee*)3/2 



Eqs. (5.23) and (5.24) are followed by the relation 



3 (l + £)(l + e*) 



1 - ee* 



-P2. 



(5.25) 



Thus, we reformulated the problem in such a way that 
the continuum begins at zero energy and there is no need 
to apply any special technique to suppress the packet 
creation energy at large distances. As a result, still re- 
maining in the framework of AMD but having redefined 
the Brink-Bloch orbitals by introducing the degrees of 
freedom describing the wave packet spreading, we can 
consider not only scattering but also resonance states in 
the relatively low energy region. 



X exp 



(RS) +e*i?V2-Ke5V2 



(5.20) 



1 - ee* 

The kinetic energy of the wave packet now takes the form 

(0,(£,R)|f|0,(£,R)) ^ 3 (l+£)(l+£*) 

((/.r(e,R)|0r(e,R)> ^4 l-ee* 

1 [(1 + £*)R-(1 + £)S]% ^^ 

"4 {l-ee*Y }(0r(e,R)|</'r(£,R)). 

Finally, the potential energy is 

(0r(£,R)| - Foexp{-(2- l)/z r2}|0,(£,R)) 
(0r(e,R)|0r(e,R)) 



(5.21) 



A 

a: 



3/2 



(1-z) A /(l-£*)R+(l-e)S^ 



(5.22) 



I 2 Az \ 1 — ££* 

where 

A = 1 - ££*, A^ ^(1 - ££*) + z)(l - £)(1 - £*). 

These expressions suffice to derive the classical equa- 
tions of dynamics and study the problem of spreading of 
the Brink-Bloch wave packet and its influence to motion 
of the packet in the field of a Gaussian potential. 

The integrals of free motion are kinetic energy (equal 
to the total energy E) 

3 (l+£)(l + £*) _ 1 [(1 + £*)R - (1 + £)S]2 

4 1 - ££* 4 



*\2 



(1 -££*) 



and the momentum 

i (1 + £*)R- (l-h£)S 



V2 



1 - ££* 



(5.23) 
P. (5.24) 



6. CONCLUSION 

Trying to describe collisions of nuclei as collisions of 
their nucleons by the classical equations of dynamics we 
face a number of problems. Let us summarize them. 

i) The Pauli principle affects the derivation of the clas- 
sical equations in the region of low values of coordinates 
and momenta where quantum corrections are expected 
to be significant. However, the corrections due to the 
Pauli principle do not suffice to describe some important 
features of wave functions by phase trajectories if the en- 
ergy of relative motion of the nucleons is low. At the 
same time, as energy increases, the classical results limit 
to the quantum ones. 

ii) The kinetic energy of classical motion of nucleons 
contains constant factors which raise the threshold en- 
ergy of break-up channels as well as the total break-up 
energy. This unfortunate phenomenon can be eliminated 
by introducing those degrees of freedom which describe 
the spreading of wave packets. The spreading affects 
the process as long as the packets move slowly and their 
width can change within the time of the nuclear reaction. 

iii) In AMD, the nucleon-nucleon interaction is taken 
into account twice. First, it forms the trajectory of each 
nucleon. Second, it is responsible for the effective cross- 
section of the nucleon-nucleon collision. A similar situa- 
tion occurs when the Vlasov equation is used for the de- 
scription of a system of particles if, along with the mean- 
field potential, a collision term is introduced. This term 
should be then expressed via the scattering cross-sections 
formed by the residual interaction. It is understood that 
the Vlasov (or Boltzmann) equation is applicable only if 
the nucleon density is low. Therefore their usage in stud- 
ies of heavy-ion collisions require additional justification. 

iv) In the ground states of even-even nuclei, the AMD 
wave function features an a-cluster structure or, if the 
number of neutrons exceeds the number of protons, some 
of the clusters are singlet pairs of neutrons. It is the 
Volkov force which leads to this conclusion. Therefore, to 
describe collisions of such nuclei, effective cross-sections 
of elastic and inelastic scattering of a- and di-neutron 
clusters should be used. Although there are experimental 
data for the a — a scattering available, the cross-sections 
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of scattering of singlet di- neutrons have to be calculated 
theoretically. The nucleons themselves appear after the 
first collision only, with some probability. Besides, inelas- 
tic scattering may produce deutron (d), triton (t) and 



helionic (h) clusters. Later they all interact with alphas 
and di- neutrons. This leads to the stochastisation, and 
the process ceases to be deterministic. 
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